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ABSTRACT 

We study the effect of a single, instantaneous starburst on the dynamical 
and chemical evolution of a gas-rich dwarf galaxy, whose potential well is 
dominated by a dark matter halo. We follow the dynamical and chemical 
evolution of the ISM by means of an improved 2-D hydrodynamical code 
coupled with detailed chemical yields originating from type II SNe, type la 
SNe and single low and intermediate mass stars (IMS). In particular we follow 
the evolution of the abundances of H, He, C, N, O, Mg, Si and Fe. We find that 
for a galaxy resembling IZwl8, a galactic wind develops as a consequence of 
the starburst and it carries out of the galaxy mostly the metal-enriched gas. In 
addition, we find that different metals are lost differentially in the sense that 
the elements produced by type la SNe are more efficiently lost than others. 
As a consequence of that we predict larger [a/Fe] ratios for the gas inside 
the galaxy than for the gas leaving the galaxy. A comparison of our predicted 
abundances of C, N, O and Si in the case of a burst occurring in a primordial 
gas shows a very good agreement with the observed abundances in IZwl8 as 
long as the burst has an age of ~ 31 Myr and IMS produce some primary 
nitrogen. However, we cannot exclude that a previous burst of star formation 
had occurred in IZwl8 especially if the preenrichment produced by the older 
burst was lower than Z = 0.01 Zq. Finally, at variance with previous studies, 
we find that most of the metals reside in the cold gas phase already after few 
Myr. This result is mainly due to the assumed low SNII heating efficiency, and 
justifies the generally adopted homogeneous and instantaneous mixing of gas 
in chemical evolution models. 

Key words: galaxies: individual: IZwl8 - hydrodynamics - ISM: abundances 
- ISM: bubbles. 



1 INTRODUCTION 

Dwarf irregular galaxies (DIG) are playing an increas- 
ingly central role in understanding galaxy evolution. 
This kind of galaxies generally has a low metallicity 
(from 0.5 Zq to 0.02 Zq), a high gas content (up to 
~ 10 times the stellar content) and their stellar pop- 
ulations appear to be mostly young. All these features 
indicate that these galaxies are poorly evolved objects. 

Many gas-rich dwarf galaxies are known to be in a 
starburst phase, or are believed to have experienced pe- 
riods of intense star formation in the recent past. These 
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galaxies are generally called blue compact dwarf (BCD) 
galaxies (Sandage & Binggeli 1984). In general, dwarf 
gas- rich galaxies, given their simple structures and small 
sizes, are excellent laboratories to investigate the feed- 
back of starbursts on the interstellar medium (ISM), 
and to study their chemical evolution. The aim is to re- 
produce the observed abundance ratios, to trace their 
recent star formation history and to discover if these 
galaxies could be the source of the intracluster gas (Gib- 
son & Matteucci 1997). 

Many authors have tried to connect late-type gas- 
rich (DIG and BCD) and early-type gas-poor dwarf 
galaxies (dwarf ellipticals and dwarf spheroidals) in 
an unified evolutionary scenario. The favourite theory 
about ISM depletion in gas-rich dwarf galaxies is based 
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on the starburst-driven mass loss. The basis of this 
model, proposed by Larson (1974) and then applied 
specifically to dwarfs by Vader (1986) and Dekel & Silk 
(1986), is that the ISM is blown out of the galaxy by 
the energetic events associated with the star formation 
(stellar winds and supernovae). The well-known correla- 
tion between mass and metallicity found for both late- 
type and early-type dwarf galaxies (Skillman, Kennicut 
& Hodge 1989) is a natural result of the increasing in- 
ability of massive galaxies to retain the heavy elements 
produced in each stellar generation. At the present time 
is not yet clear if galactic winds are really the key point 
for understanding the formation and evolution of dwarf 
galaxies (see Skillman & Bender 1995 and Skillman 1997 
for critical reviews about this point), but they certainly 
play an important role, regulating the mass, metal en- 
richment and energy balance of the ISM. 

Observational evidences in support of the presence 
of outflows have been found recently in a lot of gas- 
rich dwarf galaxies, like NGC1705 (Meurer et al. 1992), 
NGC1569 (Israel 1988), IZwl8 (Martin 1996) and many 
others. In their search for outflows in dwarf galaxies, 
Marlowe et al. (1995) pointed out that this kind of phe- 
nomena is relatively frequent in centrally star-forming 
galaxies. Again, they note a preferencial direction of 
propagation along the galaxy minor axis. In spite of 
these observational evidences, it is often difficult to esta- 
blish if the gas will leave definitively the parent galaxy. 
In order to understand the final fate of both the swept- 
up gas and the metals ejected during the starburst 
and to study possible links between early and late-type 
dwarfs, numerical simulations are needed. 

There are a lot of recent hydrodynamical simula- 
tions concerning the behaviour of the ISM and the met- 
als ejected by massive stars after a starburst. These sim- 
ulations generally agree on the fact that galactic winds 
are not so effective in removing the ISM from dwarf 
galaxies, but disagree on the final fate of the metal- 
enriched gas ejected by massive stars. Many authors 
(D'Ercole & Brighenti 1999, hereafter DB; MacLow & 
Ferrara 1999, hereafter MF; De Young & Heckman 1994; 
De Young & Gallagher 1990) have found that galac- 
tic winds are able to eject most of the metal-enriched 
gas, preserving a significant fraction of the original ISM. 
Other authors (Silich & Tenorio-Tagle 1998; Tenorio- 
Tagle 1996) have suggested that the metal-rich mate- 
rial is hardly lost from the galaxies, since it is at first 
trapped in the extended haloes and then accreted back 
on to the galaxy. 

However, all these models consider only the effects 
of stellar winds and SNII explosions on the dynamics of 
the ISM. In this paper we present models which take 
into account also the energetic contribution and the 
feedback from intermediate-mass stars and SNela, us- 
ing the most up-to-date supernova rates. The effect of 
SNIa explosions is certainly fundamental for the late dy- 
namical evolution of the ISM (up to ~ 500 Myr after 
the burst), even if their number is small. 

There is an extensive literature about the chemical 
evolution of starburst and blue compact dwarf galax- 
ies (see e.g. Matteucci & Chiosi 1983; Matteucci & Tosi 
1985; Olofsson 1995). Pilyugin (1992, 1993) and Mar- 



coni et al. (1994) suggested the idea that the spread 
in the chemical properties of these galaxies, in partic- 
ular the observed spread in He/H vs. O/H and N/O 
vs. O/H, could be due to self-pollution of Hn regions 
coupled with 'enriched' or 'differential' galactic winds. 

Therefore it is interesting to test the differential 
wind hypothesis with an hydrodynamical approach. In 
our models we are able to follow the evolution in space 
and time of the abundances of several chemical elements 
(H, He, C, N, O, Mg, Si, Fe); in particular we follow, 
with suitable tracers, the gas released by stars of dif- 
ferent initial mass. The chemical composition of each 
of these tracers is obtained by adopting the nucleosyn- 
thesis prescriptions from various authors (Woosley & 
Weaver 1995, hereafter WW; Renzini & Voli 1981, here- 
after RV; Nomoto, Thielemann & Yokoi 1984, hereafter 
NTY). 

In section 2 we describe the model and the assump- 
tions adopted in our simulations. The results are pre- 
sented in section 3 and compared with the observational 
constraints available for the BCG IZwl8. A discussion 
is presented in section 4 while some conclusions and fu- 
ture improvements of the model are discussed in section 
5. 



2 ASSUMPTIONS AND EQUATIONS 

2.1 The gravitational potential and the gas 
distribution 

It is convenient, for computational reasons, to model 
BCD galaxies. In these galaxies, in fact, the starburst 
occurs near the optical centre and the ISM structure 
is highly axisymmetric. In particular, we will focus on 
the galaxy IZwl8 which is a well-studied, very metal- 
poor BCD galaxy. IZwl8 shows very blue colors (17 — 
B — —0.88, Van Zee et al. 1998), which are indicative 
of a dominating very young stellar population, although 
one cannot exclude an underlying older one (Aloisi et 
al. 1999). Therefore, IZwl8 is an excellent candidate to 
compare with a single-burst model, although our model 
cannot reproduce the real galaxy in detail. 

Many ingredients play an important role in the dy- 
namical evolution of the ISM: the galactic structure 
(stellar component, gaseous component, dark halo), the 
energy and mass injection rate of newly formed stars 
and the size of the starburst region. 

We model the ISM of IZwl8 assuming a rotating 
gaseous component in hydrostatic isothermal (T g = 10 3 
K) equilibrium with the galactic potential and the cen- 
trifugal force. The potential well is the sum of two 
components. The first is given by a spherical, quasi- 
isothermal dark halo truncated at a distance r th , in or- 
der to obtain a finite mass: 



ph(r) = pho 



1+ - 

\r c h 



(1) 



where r = V R 2 + z 2 and r c h is the core radius of the 
dark component (we are using cylindrical coordinates) . 
According to values found in literature for the total mass 
of IZwl8 (Lequeux & Viallefond 1980; Van Zee et al. 
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1998), the halo mass is assumed to be 6.5 x 10 Mq. 
Since we do not take into account the self gravity of 
the gas, in order to reproduce the oblate distribution 
of gas inside IZwl8 (Van Zee et al. 1998), we introduce 
a fictitious 'stellar' component described by an oblate 
King stellar profile: 



LOG DENSITY 



p*{R,z) = p* 



1 + 



R 



Zc* 



(2) 



where R c * and z ci , are the core radius along the i?-axis 
and the z-axis respectively. This profile is truncated at 
the tidal radii _R t * and z t *, in order to obtain a finite 
mass M* = 6 x 10 5 M©. This structure is flattened along 
the z-axis and we assume R ci ,/z cir = Rt*/z t * = 5 and 
Rt*/R c * = 2t*/zc* = 4.29. All the structural parameters 
of our galactic model are summarized in Table 1. The 
atomic number density of the neutral ISM is defined 
as n g = 2m ^ h , where p is the ISM mass density and 
p — 7/11 is the mean mass per particle of the fully 
ionized gas, assuming a primordial abundance. 

Although this structure is rather flat, its potential 
is rounder. The gas settled in such a potential assumes 
an oblate structure resembling that of the ISM of IZwl8 
in a region R < 1 Kpc and z < 730 pc, which we call 
'galactic region'. We note however that the elongation 
is also due to the assumed rotation of the gas which is 
responsible for the flaring at large radii (see Fig. 1, up- 
per panel). Details of how to build such an equilibrium 
configuration can be found in DB. The lower panel in 
Fig. 1 shows the resulting column density of the ISM. 

We ran several models varying the gas mass and the 
burst luminosity. We discuss in detail three of them, Ml, 
M2 and M3 (see Table 2). We also describe model MC, 
similar to Ml, in which heat conduction is allowed. 



2.2 The equations 

To describe the evolution of the gas we solve the time- 
dependent, Eulerian equations of gasdynamics with 
source terms, that we write in the form: 

dp 
dt 



+ V • (pv) = ap,, 



— + V- {q v) = a p„, 



+ V • (m (g) v) = pg — (7 — l)Ve + ap*t>*, 



dm 

~dT 

de 

— + V-(ev) = -(7-l)eV-w-L + apJ eo + ^v' 



(3) 
(4) 
(5) 
,(6) 



where p, m and e are the density of mass, momentum 
and internal energy of the gas, respectively. The param- 
eter 7 = 5/3 is the ratio of the specific heats, g and 
v are the gravitational acceleration and the fluid veloc- 
ity, respectively. The source terms on the r.h.s. of equa- 
tions (3)-(6) describe the injection of total mass and 
energy in the gas due to the mass return and energy 
input from the stars. In our simulations the burst is lo- 
cated in the centre of the galaxy, therefore both energy 
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Figure 1. Upper panel: initial gas density profiles. The 
density scale is logharitmic varies linearly between -33 and 
-24. Lower panel: column density of the initial ISM seen edge- 
on (dashed line) and face-on (solid line). 



sources (SNell and SNela) and mass return are concen- 
trated inside a small central sphere of ~ 40 pc of radius. 
We treat both sources as continuous, although the SNIa 
rate is rather low (~ 1.6 Myr - , see section 2.3.1). How- 
ever, an a posteriori analysis of our results, following 
Mac Low & McCray (1988), reveals that the continuous 
energy input assumption is still valid during the SNIa 
stage. V being the volume of the burst, p* = M huTBt /V, 
where Mb urs t is the total mass of stars formed during 
the burst (see next section), v* is the circular velocity of 
these stars, and a(t) = a*(t) + asNii(£) + asNia(i) is the 
sum of specific mass return rates from stars and SNe, 
respectively (see next section), eo is the injection energy 
per unit mass due to the stellar random motions and to 
SN explosions (see next section). Finally, L = n c n p A(T) 
is the cooling rate per unit volume, where for the cool- 
ing law A(T) we follow the approximation to the equi- 
librium cooling curve given by Mathews & Bregman 
(1978). 

g 1 represents the mass density of the i element, and 
a' the specific mass return rate for the same element, 
with 1 a% ~ a - ^1- (4) represents a subsystem of N 
equations which follow the hydrodynamical evolution of 
N different ejected elements (namely H, He, C, N, O, 
Mg, Si and Fe). This enables us to calculate the abun- 
dance of the ejecta relative to the pristine ISM. 

To integrate numerically the eqs. (3)-(6) we used a 
2-D hydrocode, based on the original work of Bedogni & 
D'Ercole (1986). We adopted a non- uniform cylindrical 
axisymmetric grid whose meshes expand geometrically. 
The first zone is AT? = Az — 5 pc and the size ratio 
between adiacent zones is 1.03. 
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Table 1. Galactic parameters 



M*(M ) 


M ha i o (M ) 


R c *(pc) 


-Rt*(pc) 


»*ch(pc) 


rth(kpc) 


6 x 10 5 


6.5 x 10 8 


233 


1000 


700 


10 



Table 2. ISM parameters 



Model 


M bul . st (M ) 


M g (M ) 


n g0 (cm 3 ) 


-Eb(erg) 


Ml 


6 x 10 5 


1.7 x 10 7 


1.81 


1.5 x 10 53 


M2 


3.6 x 10 5 


1.7 x 10 7 


1.81 


1.5 x 10 53 


M3 


6 x 10 5 


4.6 x 10 6 


0.49 


4 x 10 52 



M g is the ISM mass inside the galactic region defined in the text, Mb urs t is the mass of the stars formed during the burst and 
rigo is the central atomic number density. E b is the binding energy of the gas inside the galaxy. 



2.3 The star-burst 

For the sake of simplicity we focus on a single, instan- 
taneous, starburst event located at the centre of the 
galaxy. The stars are all born at the same time but 
they die and restore material into the ISM according to 
their lifetimes. 

Here we will describe the main assumptions about 
the initial mass function (IMF), stellar lifetimes and nu- 
cleosynthesis prescriptions adopted to calculate a» and 
asN- 

2.3.1 Mass return 

In order to obtain the number AN of stars with initial 
masses in the interval AM, we adopt the Salpeter (1955) 
initial mass function (IMF) 4>{M) — j^g assumed to be 
constant in space and time: 

(P(M)dM = BM~ {1+x) AM, 



aSNIl(t) = i?SNIl(t)AM/M b urst. 



(10) 



(7) 



where x = 1.35, and B is the normalization constant 
obtained from: 

M<j>(M)AM = M burst , (8) 



/ 

JO. 3 



With Mburst = 6 x 10 5 M we get a mass of ~ 1.5 x 10 5 
M for the stars with masses larger than 2 M , in agree- 
ment with the estimate of the stellar content in IZwl8 
by Mas-Hesse & Kunth (1996). Since the stellar yields 
are calculated only for stellar masses not larger than 40 
M (WW), we adopt this value as an upper limit in eq. 
(8). Given the very low number of stars more massive 
than this limit, the chemical and dynamical evolution 
of the gas is not affected by this choice. 

We assume that all the stars of initial mass between 
8 and 40 solar masses end their lifecycle as type II su- 
pernovae. The SNII rate is defined as: 



Rsmi(t) = <f>(M)\M\, 



(9) 



where M represents the mass of the dying stars at the 
time t. The mass return rate from SNII is then given 
by: 



Here AM is the mass restored into the ISM by a star 
of initial mass M, and is defined as M — M rom , where 
Mrem is the mass of the stellar remnant. 
In terms of single elements we have: 



asNii(t) = Rsmi{t)AMi/M b 

urst j 



(11) 



where AM, is the mass restored by a star of mass M in 
the form of the specific element i. 

The specific mass return from stars with M < 8 
M is given by: 



a,(t) = <f>(M)\M\AM/M bulsU 
and: 

al(t) = <P(M)\M\AMi/M huiat , 



(12) 



(13) 



where, again, M is the mass of the dying stars at the 
time t. 

To calculate the time derivative of the mass in eqs. 
(9), (12) and (13) we adopt the stellar lifetimes given 
by Padovani & Matteucci (1993): 



t(M) = 



f 1.2M -1 - 86 + 0.003 Gyr if M > 8 M 
if M < 8 M , 

[o. 334-^/1. 79-0. 2232x(7.764-log(M)) 



\ 10 f(M) Gyr 



(14) 



where f(M) 

To obtain the quantity AM appearing in eq. (10) 
and eq. (12) we took into account the results of WW 
for massive stars (M > 10 M ) and RV for low and 
intermediate mass stars (0.8 < M/M < 8), which give 
the mass restored into the ISM by the stars at the end 
of their lifetime. For the range 8 M < M < 10 M we 
have adopted suitable interpolations between the previ- 
ous two sets of data. 

In WW the total ejected masses (processed and un- 
processed) are given for each chemical element. In gen- 
eral, however, in nucleosynthesis papers only the 'yield' 
is given, namely the fraction in mass of a given element 
i which is newly formed and ejected by a star of initial 
mass M, the quantity P»m- In this case, the ejected total 
masses are computed in the following way: 
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AMi — AMXi + MPiM , (15) 

where Xi is the original abundance of the element i in 
the star. This is the case of the yields of RV. 

From the tables of WW (which contain also the 
products of explosive nucleosynthesis) and RV we have 
derived several relations between the initial stellar mass 
and the mass restored into the ISM in the form of chem- 
ical elements for single stars of masses between 0.8 and 
40 Mo, obtained by fitting the tabulated values with 
an eighth degree polynomial. The results are shown in 
Figg. 2-4 for different initial chemical compositions and 
different mixing lenght parameters. This enables us to 
obtain the temporal behaviour of a\{t) and ag NII (t) 
for each element i. The total mass ejection rates ob- 
tained by summing over all the chemical elements are: 
asNii(i) oc i~ 0,27 and a*(t) oc i~ 1,36 (see also Ciotti et 
al. 1991). 

Finally, in analogy with eq. (10) we define the spe- 
cific mass return from the SNela as: 

aSNIa(t) = 1.4i?SNIa(t)/M bm , t , (16) 
and: 

«SNIa(*) = i?SNIa(t)AM l /M burst , (17) 

where the mass ejected by each SNIa is assumed to be 
1.4 M© (the Chandrasekhar mass). According to the sin- 
gle degenerate model (SD) , SNe la are assumed to origi- 
nate from C-0 white dwarfs in binary systems which ex- 
plode after reaching the Chandrasekhar mass as a con- 
sequence of mass transfer from a red giant companion. 
This kind of supernova explosion occurs only after the 
death of stars of initial mass less or equal than 8 M©, 
which is ~ 29 Myr after the burst. iZsNia(^) is given by 
the following formula, obtained by the best-fitting of the 
SNIa rate computed in detail numerically by the model 
of Bradamante et al. (1998), when applied to the case 
of a single starburst: 

i?SNia(t) = 4.2xlO- 9 f^tij yr- 1 , (18) 

where tg is the time expressed in Gyr. 

It is worth noting that the SNIa rate in BCG is 
practically unknown and therefore it is very difficult 
to choose the right fraction of binary systems, in the 
IMF of such galaxies, of the type required to originate 
a SNIa. The rate of eq. (18) corresponds to the rate of 
Greggio & Renzini (1983) for a starburst with a frac- 
tion of binary systems A = 0.006. This rate switches on 
somewhat more gradually than in our approximation, 
reaching a maximum after ~ 40 Myr (see Greggio & 
Renzini, fig. 1), but this difference has no consequences 
in the dynamical evolution of our models. To show this, 
we ran a model (not shown in this paper), up to ~ 80 
Myr, using the rate computed by Greggio & Renzini and 
the differences with the results shown in section 3 were 
negligible. With our assumed rate, SNela contribute by 
~ 60 per cent of the total iron production after 15 Gyr, 
in agreement with predictions for the solar neighbour- 
hood (Matteucci & Greggio 1986). 

In summary, stars in different mass ranges con- 
tribute to galactic enrichment in a different way: 



Table 3. Nucleosynthesis prescriptions 



Case 


Z 


"RV 


A 








B 





1.5 


C 


0.01 Z© 





D 


0.01 Z Q 


1.5 



Z is the abundance of the unprocessed gas. If Z = 0, we 
assume X = 0.77, Y = 0.23. The solar abundances adopted 
are taken from Anders & Grevesse (1989). 

(i) For low and intermediate stars (0.8 M Q < M < 8 
Mq) we have used the RV nucleosynthesis calculations 
for a value of the mass loss parameter r\ = 0.33 (Reimers 
1975) and the mixing lenght orv = and orv = 1.5. 
The initial chemical composition is either Z — or 
Z — 1/100 Z©. These stars mainly produce He, C, N 
and s-process elements (not considered here). In partic- 
ular, N is a 'secondary' element, namely produced from 
the original C and O present in the star at birth. There- 
fore, for zero metallicity initial chemical composition no 
N would be produced. However, there is the possibil- 
ity of producing N in a 'primary' way, namely starting 
from the C and O newly formed in the star. This is 
the case of the IMS which can produce primary N dur- 
ing the third dredge-up episode in conjunction with the 
hot-bottom burning, during the thermal-pulsing phase 
occurring when these stars are on the asymptotic giant 
branch (AGB) (case orv = 1.5 of RV). Moreover, mas- 
sive stars can also produce primary N, as suggested by 
Matteucci (1986). In the nucleosynthesis prescriptions 
of WW there is some primary N from massive stars but 
is negligible. 

(ii) For massive stars (M > 10 M Q ) we have adopted 
the case B in the WW nucleosynthesis results, focusing 
our attention on the models with Z = and Z = 1/100 
Z©. These stars are responsible for the production of the 
a-elements (O, Mg and Si) and for part of the iron. The 
stars in the mass range 8 M© < M < 10 M© produce 
mainly He and some C, N and O. 

(iii) For type la SNe we have followed the results of 
NTY adopting their model W7. In this model, every 
type la SN restores into the ISM ~ 1.4 M© of gas. Most 
of this gas is ejected in the form of Fe (~ 0.6 M©) and 
the rest is in the form of elements from C to Si. 

In Table 3 we present a brief summary of the nu- 
cleosynthesis prescriptions adopted in our models. Note 
that each of these cases can be adopted for the three dif- 
ferent hydrodynamical models. Thus, for instance, here- 
after with model M1B we intend the hydrodinamical 
model Ml with the chemical option B. 



2.3.2 Energy input 

The energy input into the ISM due to the stellar ac- 
tivity is taken into account in eq. (6) through the term 
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CASE A (Z = 0, « RV = 0) 
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Figure 2. Yields from massive and intermediate stars (data taken from WW and RV yields) together with our fits (eighth 
degree polynomial best-fittings) for the case A. M e j is the total mass of gas ejected by the star and corresponds to AM as 
defined in the text. 



eo = 3fcTb/2/i, where k is the Boltzmann constant. The 
injection temperature can be written as: 



; (a*T„ + rjnasmiTn + ??iaQSNiaTi a )/a, 



(19) 



where fcT», kTu and fcTia are the energy per unit mass 
in the ejecta of stars, SNela and SNell, respectively (see 
e.g. Loewenstein & Mathews 1987 for more details), r/n 
and ?7i a represent the efficiency with which the energy 
of the stellar explosions is transferred into the ISM for 
SNell and SNela, respectively. We assume that 10 jl erg 
of mechanical energy are produced during the explo- 
sion of both types of supernova. However, we assume 
rjn = 0.03; only 3 per cent of the energy explosion is 
available to thermalize the ISM, while the rest is ra- 
diated away. This prescription is taken from the work 
of Bradamante et al. (1998), who studied in detail the 
chemical evolution of blue compact galaxies. Actually, 
some debate is present in literature about the efficiency 
of the SNII in heating the ISM in starbursts. We thus 
run also a model with r)u = 1 which however is not suc- 
cesful in describing IZwl8, as discussed at the end of 
section 3.1.5. 



For the SNIa explosions, instead, we assume rji^ — 1 
because the SNRs expansion occur in a medium already 
heated and diluted by the previous activity of SNII. 

It is worthwhile to note that we neglected the en- 
ergetic contribution of stellar winds from massive stars, 
according to the results of Bradamante et al. (1998) who 
showed that the injected stellar wind energy is negli- 
gible relative to the SN energy in this kind of galax- 
ies and to the results of Leitherer et al. (1999). This 
is mostly a consequence of the low initial metallicity 
adopted (Z — or Z = 0.01 Zq), because the mass loss 
from stars strongly depends on their metallicity (see e.g. 
Portinari et al. 1998 and references therein). 



3 SIMULATIONS 
3.1 Dynamical results 

Our reference model is Ml. In addition, we run two 
other models, M2 and M3, which have the same poten- 
tial well (see Table 2). In model M2 the mass M bur st of 
gas turned into stars is halved in comparison with Ml. 
In this case two competing effects are expected. On one 
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CASE B (Z = 0, a RV =1.5) 
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Figure 3. Yields from massive and intermediate stars for the case B. 



hand, half of the metals are produced and the result- 
ing increase in the metallicity of the ISM is expected 
to be lower; on the other hand, the galactic wind lumi- 
nosity powered by SNe is also halved, stellar ejecta are 
expelled less effectively from the galaxy and the enrich- 
ment of the galactic ISM tends to be higher. Two com- 
peting tendencies are also present in model M3 which 
has nearly one fourth of the ISM mass in comparison 
with Ml. In this case the stellar ejecta mixes with less 
gas and the metallicity of this gas is thus expected to 
become higher than in Ml, but in this case the wind 
is favoured and less metals are retained by the galaxy, 
resulting in a lower chemical enrichment. We also show 
model MC, similar to Ml, to study the action of heat 
conduction. 

In order to discuss the hydrodynamical behaviour 
of the gas we recall briefly a few results about bub- 
ble expansion in stratified media (Koo & McKee 1992, 
and references therein) . The freely expanding wind pro- 
duced by the starburst interacts supersonically with 
the unperturbed ISM thus creating a classical bubble 
(Weaver et al. 1977) structure in which two shocks are 
present. The external one propagates through the ISM 
giving rise to an expanding cold and dense shell, while 
the inner one thermalizes the impinging wind produc- 



ing the hot, rarefied gas of the bubble interior. The 
shocked starburst wind and the shocked ISM are sep- 
arated by a contact discontinuity. The density gradient 
of the unperturbed ISM being much steeper along the 
z axis, the expansion of the outer shell occurs faster 
along this direction. A bubble powered by a constant 
wind with velocity Vy,, mass loss rate M and mechani- 
cal luminosity L w = 0.5MV^ is able to break out from 
a gaseous disc if L w > 3Lb, where the critical luminos- 
ity is Lb = 17.9poH^ g Co. Co is the sound speed of the 
unperturbed medium and po is its central density. H e g 
is the effective scale length of the ISM distribution in 
the vertical (z) direction and is defined as: 



Po Jo 



pdz. 



(20) 



In our models H e g ~ 300 pc. If the wind luminosity 
is larger than L n ~ 0.35mLb, where m = VW/Co, the 
wind blows directly out of the planar medium, at least 
in directions close to the axis. If, instead, L w < L n 
the formation of a jet is possible, in which the wind is 
shocked and then accelerated again to supersonic speeds 
through a sort of de Laval nozzle created by the shocked 
ambient medium. Kelvin-Helmholtz instabilities tend to 
distort the nozzle, and stable jets can exist only for /3 = 
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Ch/Co < 30, where Ch is the cavity sound speed (Smith 
et al. 1983). 

3.1.1 Model Ml 

For model Ml L b = 2.8 x 10 36 erg s _1 . Actually, L w 
is not constant in our simulations. However, as a repre- 
sentative value, for the wind powered by SNell we have 
t w ~ 2x 10 38 erg s" 1 (cf. Fig. 9) and m ~ 300, thus the 
bubble carved by this wind is able to break out. Note 
that L w < L n , so that a jet-like structure is expected. 
As shown in Fig. 5, a jet actually propagates with a 
shock velocity K~3x 10 6 cm s~ x after 30 Myr. This 
figure also shows that the bubble shell on the simmetry 
plane has reached the maximum allowed value (Koo & 
McKee 1992) 7? max = 0.72// off (L w /L b ) 1/6 ~ 440 pc. 

The SNII wind lasts a relatively short time (29 
Myr) , and is then replaced by a weaker SNIa wind with 
L w ~ 2 x 10 37 erg s _1 (cf. Fig. 9) and m ~ 200. The 
existing jet cannot be sustained by this wind and is in- 
hibited before breaking out. The bubble as a whole stops 
to grow, and the incoming shocked wind pushes a large 
fraction of the hot SNII ejecta against the dense and 
cold cavity walls. Thus most of these ejecta would be 
located close to the cavity edge. The thermal evolution 



of these ejecta is difficult to asses. Given the spread of 
contact discontinuities due to the numerical diffusion, 
the ejecta partially mixes with the cold wall of the cav- 
ity, so that a large fraction of these metals cools off (cf. 
Fig. 14). Actually, as discussed by DB, several physi- 
cal processes, such as thermal conduction and turbulent 
mixing, produce a similar effect. In section 3.2 we con- 
sider explicitly heat conduction, but neglect the turbu- 
lent mixing (Breitschwerdt & Kahn 1988, Kahn & Bre- 
itschwerdt 1989, Begelman & Fabian 1990, Slavin, Shull 
& Begelman 1993) which is very complex and nearly im- 
possible (and probably even meaningless out of a fully 
3D simulation) to implement into the code. 

The bubble inflated by the SNIa wind is likely to 
break out through a nozzle. We note, however, that the 
gas distribution in front of the bubble along the z di- 
rection is modified by the expansion of the outer shock 
generated by the previous SNII activity; although the 
powerful SNII wind is ceased, this shock continues to 
expand with increasing velocities (V s ~ 10 7 cm s _1 at 
t — 342 Myr) because of the steep gradient of the un- 
perturbed ISM density profile. At these shock velocities, 
the post-shock gas cools quickly and its temperature 
is T = 10 3 K (the minimum allowed in our computa- 
tions) everywhere with the exception of a 'rim' behind 
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the shock, where T ~ 3 x 10 K. The density gradient 
of the upwind gas experienced by the SNIa bubble is 
shallower, and the break out is contrasted. This can be 
seen in Fig. 5, where the hot bubble is shown to grow 
very little up to t ~ 300 Myr. 

The gas in the bubble radiates inefhcently because 
of its low density (n~3x 10~ 4 cm~ 3 ) and its tem- 
perature is ~ 2 x 10 6 K. Shears are present at the con- 
tact surface between hot and cold gas, which is thus 
Kelvin-Helmholtz unstable. For this reason the cavity 
is irregularly shaped at this stage. 

Subsequently, since the surrounding cold gas is in 
expansion, the hot gas is finally able to break out carv- 
ing a long tunnel. This tunnel has the de Laval nozzle 
structure, with the transverse section increasing with z. 
We note the presence of Kelvin-Helmholtz instabilities 
at the wall of the nozzle. This is due to the fact that 
/3 ~ 30, so the nozzle is only marginally stable. 

The shocked wind is accelerated again to supersonic 
speeds through this nozzle (velocities V ~ 4 x 10 7 cm 
s _1 and mach numbers M ~ 8). When the jet is well 
developed (t = 342 Myr in Fig. 5b), the minimum radius 
of the nozzle is R n ~ 100 pc. 

The acceleration of the cold shell in front of the 
jet causes it to be disrupted by the Rayleigh- Taylor in- 
stabilities, and the hot gas leaks out. Contrary to the 
previous works where the SNIa activity was not consid- 
ered, at these late times the central galactic region is 
not yet replenished by the cold surrounding gas. Tak- 
ing into account equation (18), we stress that this will 
happen after ~ 2 Gyr, when L w = 2 Lb- 

In Fig. 9 we have plotted mass, energy and lumi- 
nosity budget inside the galaxy. From the central panel 
of this figure we note that the energy of SNe never be- 
comes larger than the binding energy, although some 
gas is definitively lost from the galaxy, as it is appar- 
ent from the numerical simulation. This indicates that 
ballistic arguments cannot be adopted properly to cal- 
culate ejection efficiencies. In fact, an element of fluid 
can acquire energy at the expense of the rest of the gas 
through opportune pressure gradients, thus increasing 
its velocity beyond the escape velocity. 

The thermal energy shows in particular two drops 
at t ~ 30 Myr and at t ~ 160 Myr. The second drop 
reflects a decrease in the hot gas content, while the 
first one cannot be associated at any particular hot gas 
loss. In fact this drop coincides with the discontinuity 
SNII/SNIa, when the specific energy injection falls by a 
factor ~ 10. Thus the thermal content of the bubble de- 
creases via radiative cooling, although its temperature 
remains larger than 2 x 10 4 K, which is the threshold 
adopted to define hot regions in the upper panel of Fig. 
9. The fall-off at t ~ 160 Myr is instead due to the 
presence of large eddies which move part of the hot gas 
outside the galaxy. We finally point out that at the be- 
ginning of the SNII activity the X-ray emission is absent 
(see lower panel of Fig. 9) because the energy injection is 
not able to rise the cavity temperature over T — 7 x 10 5 
K, which is the threshold adopted to define the X-ray 
emitting gas. 



3.1.2 Model M 2 

In this model, Lb is the same as in Ml, but L w is a 
factor of 0.6 lower. The dynamical evolution is rather 
similar to that of the previous model. Obviously the 
bubble is smaller at the end of the SNII activity. Quite 
surprisingly, however, the nozzle carved by the SNela 
breaks out earlier than in Ml. Note that in this case 
L w /Lb < 10, and the nozzle is stable and well shaped. 
Actually, this condition is equivalent to the condition 
j3 < 30 expressed above (Koo & McKee 1992) and, for 
this model, we find /3 ~ 25. Kelvin-Helmholtz instabili- 
ties are present, but they are carried away by the flow 
before they can grow significantly. Less energy is dissi- 
pated by the turbulence and is more easily channelled 
through the nozzle. 



3.1.3 Model M 3 

Because of the lower ISM density, we have Lb = 7.5 x 
10 35 erg s _1 for this model, so that L w > L n during the 
SNII stage. In this case the galactic wind breaks out 
rather vigorously, and the evolution of the gas is similar 
to that found in other theoretical works (Suchkov et al. 
1994, De Young & Heckman 1994, MF, DB). A promi- 
nent lobe is formed, which is similar to that described 
by Mac Low, McCray & Norman (1989) in their fig. 2. 
The shell accelerates and becomes Rayleigh- Taylor un- 
stable, expelling the hot interior. The breakup occurs at 
the polar cap, where the ISM has the lower density and 
pressure. The hot gas which blows out of the bubble pro- 
duces the jet-like structure visible in Fig. 7a. Note that 
the outer shock initially is rather slow (V B ~ 5 x 10 6 cm 
s _1 ), and then accelerates somewhat up to V s = 1.6x 10 7 
cm s _1 . Thus, also in this model the shocked gas never 
reaches high temperatures and most of the lobe volume 
is cold. 

Fig. 7b shows the late dynamical evolution of the 
gas. The gas flow expands along a conical configura- 
tion, inside a solid angle which remains constant during 
all the simulation. The ISM outside this funnel remains 
substantially unperturbed. Actually, the aperture of the 
cone in our model is evidently dictated by the assumed 
structure of the ISM and could be not realistic. How- 
ever, given that no falling back or fountain is expected 
by the expelled material which is kept in expansion by 
the SNela, the final chemical characteristics of the gas 
inside the galaxy are not affected by the exact structure 
of the 'chimney'. 

Concerning the distribution along the z direction, 
the lobe of outflowing gas can be grossly divided in 
three regions. The most external, far from the galaxy, is 
bounded by the outer shock and hosts mainly shocked, 
low-density external medium. The inner region is filled 
with gas ejected by SNela and low-mass stars. Between 
these two regions there is the gas of the 'first' bubble 
(where the ejecta of SNII are present), which is quickly 
cooled to low temperatures. At t ~ 150 Myr, the shell of 
the small bubble breaks, like the first one, and the hot 
gas flows forward into the lobe rising its temperature up 
to 10 6 K. 

Just before breaking out, the superbubble reaches 
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Figure 5. Density contours and velocity fields for the model Ml at different epochs. The density scale (logarithmic) is given in 
the strip on top of the figure. In order to avoid confusion, we draw only velocities with values greater than 1/10 of the maximum 
value. This is true also for Figg. 6, 7 and 8. 
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Figure 6. As Fig. 5, but for model M2. 
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Figure 7. As Fig. 5, but for the model M3. 
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Figure 8. As Fig. 5, but for the model MC. 
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the edge of the galaxy even along the i?-axis, pushing 
out all the unprocessed gas present, and almost all the 
galaxy is covered by the hot cavity. As the breakout 
occurs, the bubble shell shrinks slightly and part of it 
comes back into the galactic region producing the rise 
of the mass of hydrogen and the other elements (cf. Fig. 
11). This happens in coincidence with the pressure de- 
crease in the SNIa bubble following the rupture of the 
unstable shell, as discussed above. 

3.1.4 Model MC 

In Fig. 8 we show model MC, identical to model Ml 
but with the heat conduction activated. In order to take 
into account the thermal conduction, we solve the heat 
transport term through the Crank-Nicholson method 
which is unconditionally stable and second order accu- 
rate (see DB for more details). In this model the cavity 
is less extended than in Ml because of the increased 
radiative losses due to the evaporation front. During 
the SNII stage the bubble never extends beyond H e ff. 
Thus the "nose" present in Ml does not develop, and 
the bubble has a more round aspect (Fig. 8a). Thermal 
conduction smooths temperature inhomogeneities and 
the cavity is more regularly shaped also at later times 
(see Fig. 8b). Less energy is dissipated through eddies, 
and the final break out is slightly anticipated compared 
to model Ml. The resulting outflow is stable and well- 
shaped. The fraction of cold ejecta does not change sub- 
stantially compared to model Ml (see section 3.2 for a 
discussion about this point). 

3.1.5 ISM ejection efficiencies 

It is useful to define the efficiency of gas removal from 
the galaxy. This efficiency cannot be unambigously de- 
fined as in the case of ballistic motions because of dis- 
sipative effects which may play a very important role 
(DB). For this reason we simply define the efficiency 
/ism by dividing the mass of the gas which has left the 
galaxy by the total mass of the ISM. In a similar way we 
calculate the efficiency in the ejection of material from 
SNII (/snii), intermediate-mass stars (/ims) and SNIa 

(/SNIa)- 

The upper panels of Fig. 11 shows the masses of 
the different elements removed from the galaxy as func- 
tions of time. The relative proportions between masses 
of different elements is essentially that expected by a 
Salpeter IMF and is not substantially affected by se- 
lective dynamical losses. For the model M1B, we point 
out that the mass of metals lost from the galaxy de- 
clines after reaching a maximum at t ~ 200 Myr. This 
maximum is mirrored by a minimum in the masses of 
metals inside the galaxy. This behaviour is due to the 
fact that, at this time, the large hot blob visible in Fig. 
5b (second panel) extends over the galactic edge (R di- 
rection) thus inducing a large loss of ISM (mostly H 
and He). We calculate the efficiency at t ~ 200 Myr and 
at the end of the simulation (t ~ 375 Myr), obtain- 
ing (/ism, /snii, /sNia,/iMs)=(0.43,0.38,0.25,0.38) and 
(/ism, /snii, /sma, /ims)=(0.07,0.17,0.20,0.26), respec- 
tively. Thus, at t ~ 200 Myr the products of the SNII 




- Total mass 

Feedback from massive stars 
■ Feedback from intermediate stars 

- Feedback from SNIa 
Mass of hot regions 



■ SNIa energy 

■ SNII energy 

■ Binding energy 
Kinetic energy 

■ Thermal energy 



SNIa luminosity 
■ SNII luminosity 

At 



7 7.5 8 8.5 

log I (yr) 

Figure 9. Energy, luminosity and mass budget inside the 
galaxy for model Ml. Hot regions are defined as the regions 
where T > 2 X 10 4 K. Lj indicates the emission from gas 
with T > 7 X 10 5 K" (emission in the X-ray band), while Ltot 
is the total emission. 



and IMS have been ejected more easily than the prod- 
ucts of SNela. This is of course due to the fact that SNIa 
material is located in a region closer to the galactic cen- 
tre. After the break up all the efficiencies decrease, and 
in particular /snii shows the greatest reduction. In fact, 
the SNII ejecta inside the galaxy are 'incorporated' into 
the cold, dense shell of the cavity and do not experi- 
ence any substantial further dynamical evolution (see 
discussion in section 3.1.1 and Fig. 5b); at the break 
out the bubble shrinks and its walls recede entirely in 
the galaxy, increasing its content of SNII ejecta. Con- 
trary to the expectations, the higher efficiency is given 
by /ims instead of /sNia- This is due to strongly un- 
steady behaviours of the nozzle wall. 

In model M2, where a nearly steady flow is ob- 
tained, we actually have /snii < /ims < /sNia. 
At t ~ 300 Myr it is (/ism, /snii, /sNia, /ims) 
= (0,0.06,0.32,0.12) in this model. For this model the dif- 
ference in efficiencies between the total gas and metals 
is particulary striking and indicates that the differential 
galactic wind assumption, adopted in several one-zone 
chemical models, is a natural outcome in this scenario. 
In particular we note that at late times the galaxy is 
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almost completely replenished by gas. In fact, as ap- 
parent in Fig. 6b, the nozzle has a rather small section 
(R n ~ 85 pc) and the volume of the cavity is negligible 
in comparison with the galactic volume. The more reg- 
ular hydrodynamical behaviour reflects also in the more 
regular temporal trend of the ejected masses. Note that, 
as expected, in M2 the galactic wind starts later in com- 
parison with Ml. 

Model M3 predicts of course the maximum amount 
of metals lost and is also the first in which the break out 
occurs. The striking minimum in the metal contents of 
the galaxy occurring at t ~ 130 Myr is due to the dy- 
namical behaviour of the buble shell, as discussed above. 
The efficiencies for the model M3 at the end of the 
simulation (t ~ 470 Myr) are (/ism, /snii, /sNia, /ims) 
= (0.77,0.85,0.97,0.87). 

Finally, we mention that, as outlined in section 
2.3.2, we run a model (not shown here) similar to Ml, 
but with r/n = 1. At the end of the SNII stage, the 
galaxy is almost devoided of gas (a part for the tenuous 
galactic wind) . The time-scale for the galactic replenish- 
ment (the SNIa wind cannot preserve an empty region 
so large) is at least R t ,/Co ~ 200 Myr. Actually, due 
to the retarding effect of the centrifugal force, in our 
simulation most of the galaxy is replenished after ~ 450 
Myr. The age of the (last) burst occurred in IZwl8 is 
estimated to be < 27 Myr (see below), thus its actual 
content of gas rules out the possibility of an high r/u. 
Note that our assumption that all energy injection oc- 
curs in the central region leads to the most effective gas 
removal for a given luminosity (Strickland & Stevens 
1999). Models with a more realistic burst diffuse across 
the galaxy will be presented in a next paper. 

3.2 Instantaneous versus delayed mixing 

In the previous models a large fraction of the stellar 
ejecta cools quite soon and a rapid mixing is expected 
given the relatively short diffusion time at these temper- 
atures (for a comparison of the diffusion times in the dif- 
ferent ISM phases, see for example Tenorio-Tagle 1996). 
This is an important point in view of the confrontation 
of the results of our models with the observable abun- 
dances of IZwl8 (cf. section 4.2), and deserves some 
discussion. 

Rieschick & Hensler (2000), for instance, presented 
a chemodynamical model of the ISM of a dwarf galaxy 
in which the metal enrichment undergoes a cycle last- 
ing almost 1 Gyr. This model is based on the scenario 
described in Tenorio-Tagle (1996). In this scenario the 
break out is inhibited and the SNII ejecta (SNela are 
not considered), mixed with hot evaporated ISM, are 
located inside a large cavity which extends above (and 
below) the galactic disc. Typical parameters of this cav- 
ity are the linear size d > 1 kpc, the density n c = 10~ 2 
cm -3 and the temperature T c — 10 6 K. After the last 
SN explosion, strong radiative losses occur. However, 
given the density and temperature fluctuations in the 
hot medium, cooling acts in a differential way. This leads 
to condensation of the metal-rich gas into small molecu- 
lar droplets (_R C ~ 0.1 pc, M c ~ 1 M ) able to fall back 
and settle on to the disc of the galaxy. With the next ex- 



ploding stellar generation, the droplets are dissociated 
and disrupted, and their gas is eventually mixed in the 
HII regions. 

Here we point out some caveat concerning this sce- 
nario which should be taken into account. 

Thermal conduction - Thermal conduction, if not 
impeded by magnetic fields and/or plasma instabilities, 
introduce, together with radiative losses, the character- 
istic Field length (Begelman & McKee 1990, Lin & Mur- 
ray 2000) 



3k(T)T \ 
A| - 1 i. 2 /3<A(T) J 



1/2 



where k(T) = 6 x 10~ 7 T 25 erg cm" 1 KT 1 is the clas- 
sical thermal conductivity (Spitzer 1956) and n the 
gas density. The cooling rate scales linearly with the 
metal content £ of the gas (C — 1 f° r solar abundance). 
The parameter /3 takes into account the possibility that 
the cooling gas may be out of ionization equilibrium; 
Borkowski, Balbus & Fristrom (1990) have shown that 
13 may be as high as 10 through conductive fronts. 

Clouds undergo evaporation unmodified by radia- 
tive losses if they are sufficiently small, R c < Ap. As- 
suming T c = 10 6 K and n c = 10~ 2 cm" 3 , Tenorio-Tagle 
obtains R c < 10 pc for the radius of the overdense zones 
at the beginning of their implosion. Adopting f3 = £ = 1 
and A = 1.6 x 10" 1S T" 7 ergs cm 3 s" 1 for 10 5 K < T < 
10 7 - 5 K (Mac Low & McCray 1988), we obtain A F ~ 138 
pc. Even assuming /3 = 10, Ap ~ 44 pc remains larger 
than R c . Thus the rate of conductive heat input exceeds 
that of the radiative losses and the cloud collapse is in- 
hibited. This result derives from the general property of 
the evaporation to stabilize thermal instabilities (Begel- 
man & McKee 1990). Actually, Tenorio-Tagle obtained 
i? c ~ 10 pc as a lower limit for the radius value of clouds 
able to implode. As a consequence, perturbations with a 
size larger than A F can actually grow. A non negligible 
fraction of gas mass would condense only for a rather 
flat size spectrum of the fluctuations. 

We point out that, even if droplets actually form, 
they evaporate in a time r e = M c /M. Droplets have 
rather small final radii of order 0.1 pc (Tenorio-Tagle 
1996). In this case they undergo the satured evaporation 
M = 1.22 x 10- 14 T 2 - 5 R c o-- 5/8 g s"\ where o = 3 x 
10 18 (T/1.54 x 10 7 ) 2 /(ni? c ) is the saturation parameter 
(Cowie, McKee & Ostriker 1981). For M c ~ 1 M a we 
have r ~ 40 Myr, comparable to the dynamical time 
Td (see below) . Thus all the droplet material, or at least 
a large fraction of it, returns into the hot phase before 
reaching the galactic plane. 

Drag disruption - Suppose that thermal conduction 
is impeded and that droplets with final R c ~ 0.1 pc 
actually form and fall toward the galactic plane. During 
their descent the droplets experience a drag reaching a 
terminal speed Vt ~ (xRc/d)°' 5 V c , where \ = 10 5 is the 
ratio of the droplet density to the hot gas density, d ~ 1 
kpc is the droplet distance to the galactic plane, and V c 
is the circular velocity of the halo potential. The motion 
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of the droplet relative to the hot gas leads to mass loss 
through Kelvin-Helmholtz instability. For wavelengths 
A ~ R c the stripping time-scale is r s ~ RcX°' 5 /Vt — 
T d {R c /d) B (cf. Lin & Murray 2000), where r d = d/V c is 
the dynamical time. The droplets are thus disintegrated 
before they settle to the galactic plane, returning to the 
hot diluted phase. Larger droplets may not be able to 
attain their terminal velocity, but even in this case we 
have r s /r d ~ X°' 5 Rc/d < 1. 

As an aside, we note that, if a large fraction of the 
hot gas becomes locked into droplets, the pressure of the 
remaining diluted phase reduces and the cavity shrinks. 
Whether the bubble deflates slowly or suddenly and pro- 
ducing turbulence depends on the droplet formation ef- 
ficiency. Thus, the scenario depicted by Tenorio-Tagle 
of a nearly steady hot cavity with size of > 1 kpc wait- 
ing for the onset of radiative cooling could be incorrect. 
Part of the droplets could be overrun by the edge of the 
imploding bubble, undergoing an even faster stripping. 

Let us consider now the results shown in the present 
paper. Contrary to the scenario sketched above, the 
ejecta cool rapidly without leaving the galaxy (until the 
break-out, which occurs at late times) and without un- 
dertaking a long journey before mixing with the ISM. 
How much these results are reliable? Mac Low & Mc- 
Cray (1988) showed that a conductive bubble expanding 
in an uniform medium becomes radiative (i.e. radiates 
an energy comparable to the thermal energy content of 
the shocked wind) after a time: 

t R ~ 16(/3C)- 35/22 ^ n n- 8/11 Myr, 

when the cavity radius is: 

Rn ~ 350(P(r 27/22 LTn- 7/11 pc. 

For t > £r the bubble goes out of the energy conserving 
regime, although a fully momentum conserving regime 
is never attained. Considering model Ml, we assume 
ft = £ = 1, L38 = 2 and n = 1.8 cm~ 3 , and we obtain 
*r = 12.6 Myr and Rn — 310 pc. Thus, in the case of an 
uniform unperturbed medium the bubble interior would 
cool quite early, when it is still well inside the galaxy. 

Mac Low & McCray also considered the expansion 
in a stratified medium. Equating the radius of a spher- 
ical bubble to approximatively one scale height H they 
define the dynamical time 

t D ~ ^(pAM 173 . 

Then the ratio of cooling to dynamical time-scales is: 

g = 8.22^35/33^0/33 ^_A_y j/3 (/3cr 3 5/22i (21) 

(note that the numerical coefficient in this expression 
slightly differs from that obtained by Mac Low & Mc- 
Cray). For model Ml we obtain tn/to ~ 1-08 during the 
SNII stage. Thus a non negligible fraction of the wind 
luminosity is radiated away (see also Fig. 9), and the 
break out does not occur. For M3, where n ~ 0.5, we 
have tn/to ~ 4.2, and the situation, in principle, is less 
clear-cut (see, however, below in this section). 

Although models Ml, M2, M3 do not take explicitly 



into account heat conduction, yet they obtain results ac- 
cording to the above scenario. In fact, as pointed out in 
section 3.1.1, numerical diffusion simulates thermal con- 
duction originating spurious radiative losses which oth- 
erwise would be absent. Of course, this spurious cooling 
does not reproduce quantitatively the same amount of 
radiation lost through a real heat conduction front, and 
the fraction of cold ejecta obtained in our models could 
be larger than the correct one. Some algorithms may be 
conveniently adopted to reduce this effect. Consistent 
advection (Stone & Norman 1992) is implemented in 
our code and helps in reducing somewhat the diffusion, 
making it consistent for the all advecteded quantities 
(mass, momentum, energy). We also made tests modif- 
ing the cooling algorithm in presence of unresolved con- 
tacts, following Stone & Norman (1993). Although the 
fraction of cold ejecta reduces of 15 per cent, most of 
the metals (~ 80 per cent) remains cold. We point out, 
however, that in presence of an unresolved conduction 
front, the above algorithm may lead to an excessive re- 
duction of the radiative losses (see below). 

In any case the intrinsic diffusion of the code may be 
alleviated but not eliminated by algorithms as those de- 
scribed. In principle, more realistic models can be done 
explicitly adding the physical terms which produce dif- 
fusion. For this reason we also ran model MC, where 
the heat transfer is included. In this model the amount 
of cold metals does not change appreciably. However, 
although model MC is useful to understand the stabi- 
lizing effect of heat conduction on a turbulent flow, it 
turns out to be inadequate to obtain the correct cooling 
rate at the conduction front. Consider the temperature 
profile of a "standard" bubble Tb(l — r/rb) 2 ^ 5 (Weaver 
et al. 1977), where Tb = 10 6 K is the central tempera- 
ture, and r b = 300 pc is the bubble radius (cf Fig. 7). 
The cooling curve maximum occurs at T ~ 2 x 10 5 K. 
This temperature is found at r/r b = 0.98, i.e. at a dis- 
tance of 6 pc to the cold shell. At a distance of 300 pc 
the mesh size is ~ 15 pc, and the conduction front is 
not resolved properly. Thus, we also ran models MCH 
and MCHH (not shown here) with heat conduction and 
with an uniform grid with mesh size of 2 and 1 pc re- 
spectively. These models were computed only up to the 
end of the SNII activity because of the large number of 
grid points involved. 

The four panels in Fig. 10 show the profiles along 
the galactic plane (z = 0) of several quantities for mod- 
els Ml, MC, MCH and MCHH at nearly 30 Myr. As ex- 
pected, the resolution of the temperature profile is not 
improved in MC and the temperature jump remains un- 
resolved. In MCH, instead, this jump extends over 2-3 
meshes and in MCHH over 4-5, as expected. The frac- 
tion of cold ejecta is 0.95, 0.95, 0.93, 0.92 for Ml, MC, 
MCH, MCHH respectively. Although the greater accu- 
racy, the fraction of cold ejecta in MCH and MCHH is 
only few percent lower than in Ml. This occurs because 
the bubble is "genuinely" radiative, and an high spatial 
resolution may retard a little the cooling, but cannot 
avoid it. 

We stress once more (cf. section 3.1.1) that it is 
very difficult to give the correct description of the con- 
tact discontinuity at the outer edge of the hot cavity, 
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for the presence of complex hydrodynamical phenom- 
ena occurring there. These phenomena tend to produce 
a finite thickness of the contact, giving rise to a substan- 
tial cooling otherwise absent. The correct evaluation of 
such a cooling is however very difficult to assess. Even 
restricting ourselves on the simple case of heat conduc- 
tion, the possible lack of ionizing equilibrium and the 
numerical diffusion both increase radiative losses. The 
first effect is physical and could be described solving 
step by step the time dependent set of ionization equa- 
tion. The second is a spurious result due to numerical 
diffusion of the code which must be reduced as much 
as possible. We believe that our convergence test (cf. 
fig. 10), as well other ID tests not shown here, indicates 
that our results are significative. We are aware of the 
fact that more refined simulations may produce some- 
what different values of the cooled mass of the ejecta, 
but we think that the large fraction of resulting cooled 
ejecta is a genuine result. 

As pointed out above, equation (21) gives an am- 
biguous prevision about the behaviour of M3. We there- 
fore ran also this model on a uniform grid with 2 pc 
resolution (up to t — 30 Myr) adding heat conduction. 
The overall dynamics of the superbubble remains the 
same, and the fraction of cold ejecta turns out to be 
6 per cent lower than in the low resolution model. We 
thus conclude that the results obtained by our models 
are reliable. 

At the end of this section we mention the effects 
expected in the case of a non-uniform initial ISM. Ac- 
tually, gas clouds are embedded in the pervasive diffuse 
gas of the galaxy. However, a correct numerical treat- 
ment of the interaction between these clouds and the 
ambient gas introduces enormous complications in the 
description of the involved physics and needs full 3D 
computations with an huge number of meshes. We here 
just make some simple considerations following McKee, 
Van Buren & Lazareff (1984). These authors describe 
the behaviour of a bubble generated by an O star and 
expanding in a cloudy medium. Because of the flux of 
ionizing photons emitted by the star, the nearby clouds 
undergo photoevaporation and accelerate away through 
the rocket effect. Essentially no cloud survives up to 
a radius Rh, and the gas density inside this radius in- 
creases to a value 0.5n m , i.e. half of the mean density the 
cloud gas would have if it were homogenized. The wind 
bubble evolution depends on the value of L, = L w /Lst, 
where L S t = 1.26 x 10 36 (S , ! 9 /n m ) 1/3 and S 4 g is the rate 
at which the star emits ionizing photons in units of 
10 49 s _1 . For weak winds (L* <C 1) the bubble radius 
is smaller than and it evolves "normally" (Weaver 
et al. 1977). For moderate winds (L» ~ 1) the bubble 
expands to the edge of the cloud distribution because 
photo-evaporated gas induces radiative losses reducing 
the pressure. Finally, for L* ^> 1 the bubble rapidly en- 
gulfs a number of clouds and radiates away most of its 
internal energy. This scenario cannot be directly applied 
to a star burst as a whole. In fact clouds are also present 
inside the star formation region, partially screening the 
flux of ionizing photons escaping from this region. How- 
ever, even assuming that this effect is negligible, in 
our model L* ~ 3. We in fact computed S49 running 



remotly at the site www.stsci.edu/science/starburst99 
(Leitherer et al. 1999) a burst model tailored on that 
assumed here. The number of UV photons produced by 
massive stars remains nearly constant (S49 = 584 s _1 ) 
up to 5 Myr, and then drops as t~ 4 . For n m = 1.8 cm~ 3 
close to the galactic disc, we have Lst — 7 x 10 37 erg 
s _1 , lower than L w . Of course, along the z direction R h 
could move much further, but the ionizing flux starts to 
decline rather soon, and the cloud distribution remains 
nearly unaffected. Thus the bubble is expected to be- 
come radiative sooner than in the case of a smoothly 
distributed ISM. 

In conclusion, in the scenario of our models, most 
of the metals actually cools off in a few Myr. It is worth 
noting that this result depends essentially on the as- 
sumption of a low heating efficiency of SNIIs. In the 
model similar to Ml, but with r/u = 1 (sect. 2.3.2), the 
break out quickly occurs and the metals have no time 
to cool. In fact, following Eq. (21), in this case the bub- 
ble results to be adiabatic and not radiative. Thus the 
SN efficiency value is crucial, and a future paper will be 
devoted to it (D'Ercole & Melioli, in preparation). 

3.3 Chemical results 

In Fig. 10 is shown the evolution of the element abun- 
dances for the models M1B, M2B and M3B. In this fig- 
ure, the evolution of the masses and the abundances 
in the form of the various elements is shown. It is in- 
teresting to note that the mass of the lost metals for 
the model M3B is larger than that retained from the 
galaxy, whereas it is the contrary for models M1B and 
M2B. This means that the initial conditions, namely the 
assumed burst luminosity and gas mass in the galaxy, 
are playing a very crucial role in the development and 
evolution of the galactic wind. The abundances are cal- 
culated as [Z/H] , where Z indicates the abundance of 
the following elements, O, C, N, Mg, Si and Fe, relative 
to the solar abundances of Anders & Grevesse (1989), 
and as 12+log(Z/H), which is the notation normally 
used for the abundances in extragalactic Hn regions. 
These abundances are derived in the following way: for 
the galactic abundances we have averaged in the pre- 
viously defined galactic region (approximatively an el- 
lipsoid with major semi-axis of 1 Kpc and minor semi- 
axis of 730 pc) , whereas for the abundances leaving the 
galaxy the integral is made over the rest of the grid. 

The a-elements show a very similar evolution and 
this is due mostly to their common origin. The a- 
elements are, in fact, mainly produced by massive stars 
(see Figg. 2-4), and thus their abundances inside the 
galaxy grow in the first 6-7 Myr - the time interval 
where most of massive stars die - then show a slowly 
decreasing trend, with however a little maximum around 
100 Myr. This behaviour is related to the dynamics of 
gas flows described before. 

The evolution of iron is not too different with re- 
spect to the evolution of the a-elements although this 
element is substantially produced by SNela. This is due 
to the fact that at the end of our simulations the iron 
produced by SNIa is only ~ 30 per cent of the total, 
because the bulk of SNela appears at later times. The 
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R (Kpc) 

Figure 10. Density (heavy solid lines), ejecta (dashed lines) and temperature (light solid lines) profiles for models Ml (first 
panel starting from top), MC (second panel), MCH (third panel) and MCHH (bottom panel) near the conduction front, along 
the _R direction, after ~ 30 Myr. The diamonds superimposed to the temperature profiles, indicate the mesh points of the grid 
for each numerical simulation. 



evolution of N shows a sharp increase at around 29 Myr 
which corresponds to the lifetime of a 8 M Q star, which 
is the first star producing a substantial amount of N 
(of primary origin). For times less than 29 Myr the N 
production is negligible. 



4 DISCUSSION 

We want to compare now our models with observational 
data found in literature for the galaxy IZwl8. 

4.1 Morphology and dynamics 
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Figure 11. Time evolution of several quantities for the models M1B, M2B and M3B (first, second and third column, re- 
spectively). The upper pannels show the evolution of the mass of various elements outside the galaxy, while the masses inside 
the galaxy are shown by the second raw of pannels. The third and the last raws of panels illustrate the behaviour of the ISM 
abundances (relative to the sun and by number, respectively). 



To compare our dynamical results we should first 
summarize the structural and dynamical properties of 
IZwl8. IZwl8 has a 'peanut-shaped' main body, con- 
sisting of two starbursting regions (Dufour et al. 1996). 
There are also two H n regions (also called NW and SE) , 
associated with the main body, but shifted ~ 1 arcsec 
east of the brightest continuum emission (Martin 1996). 
The Ha emission is bipolar-shaped along a direction or- 
thogonal to the main body, and show clear evidences of 



shell structures. In fact a prominent shell stretches 15 
arcsec (720 pc) north-northeast from the northwest H n 
region and bright Ha emission extends symmetrically 
south-southwest from the NW region (Martin 1996). 
Moreover a partial shell of 3.6 arcsec of diameter (173 
pc) protrudes from the north-west side. The Hi veloc- 
ity field (Van Zee et al. 1998; Viallefond et al. 1987) 
shows a significant velocity gradient along minor axis, 
suggesting a flow in this direction (Meurer 1991). 
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Table 4. Predicted abundances in the galactic region after 31 Myr for models Ml, M2, M3. We also show the results for model 
MC (with nuclcosynthctic prescriptions B), in order to emphasize the similarities with model M1B. Only models Ml and MC 
should be compared with IZwl8. These values are compared with some abundances found in literature for IZwl8. 
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We first note that the distance of 720 pc between 
the shells and the NW Hn region is quite compatible 
with models Ml, M2, M3 (in these models, after 31 Myr, 
the bubble has covered a distance of ~ 600-800 pc along 
the z-axis). Martin, starting from geometrical consider- 
ations, found a shell speed of 35-60 Km s _1 , and in our 
model Ml, after 31 Myr, the velocity of the outer shock 
along the «-axis is approximatively 30 Km s _1 , in agree- 
ment with the observations. We note instead that the 
model MC does not develop an outflow along the z di- 
rection (see Fig. 8a), and this is due to the fact that the 
bubble never extends beyond H c f{, as explained in sec- 
tion 3.1.4. Actually, the uncertainties on the real value 
of H c h in IZwl8 are significant and a small reduction of 
H e g could lead to the formation of an outflow also for 
model MC. 

These comparisons between observation and theory 
depend strongly on the adopted distance of IZwl8. In a 
recent work, Izotov et al. (1999) found that the distance 
of IZwl8 should be at least 20 Mpc, twice the distance 
generally adopted for this galaxy. With this new dis- 
tance, Color-Magnitude Diagram (CMD) studies give 
an age, derived from the main sequence turn-off, of 5 
Myr for the main body. They suppose that the star for- 
mation in the main body has started ~ 20 Myr ago in 
the NW edge, propagating then toward the SE direction 
and then triggering the main starburst ~ 5 Myr ago. 
This estimate is also consistent with the stellar popula- 
tion analysis of Hunter & Thronson (1995). 

Here we have assumed a coeval stellar population, 
so we cannot correctly verify the hypothesis of Izotov 
et al., but in our models a slightly pre-enriched burst 
with an age of 5 Myr cannot account for the observed 



abundances and the nearly flat metallicity gradient ob- 
served in IZwl8 (see next section). Moreover, with the 
new estimated distance the above-mentioned shell struc- 
tures double their dimensions and it is hard to reproduce 
these morphological features with a burst 5 Myr old. 



4.2 Chemical abundances 

In Table 4 the abundance ratios log(C/0), log(N/0), 
log(Si/0) as well as 12+log(0/H) predicted by our 
models are reported. At the beginning of section 3 we 
described several factors which affect the metal content 
of the ISM. From an inspection of Fig. 11 and Table 
4 we conclude that a reduction in the burst luminosity 
produces a reduction in the total abundances, while a 
decrease in the ISM mass leads to an increase of the 
metallicity. 

In Table 4 the observed abundances of IZwl8 
are also reported. However, only the model Ml (cases 
A,B,C,D) should be compared with IZwl8, since the 
total mass and the gas mass of this model have been 
chosen to match this galaxy. Before to continue the dis- 
cussion on our results, we point out that the abundances 
shown in Table 4 refer to the whole gas into the galaxy, 
while the comparison with the data should be valid only 
for the cold (T < 2 x 10 4 K) phase. In fact, chemical 
composition of stellar winds and supernovae ejecta are 
mainly measured through the relative intensities of vi- 
sual [On], [Oni], [Sn] and [Nil] forbidden lines com- 
pared to H and He recombination lines, but this ap- 
proach is sensitive only to emission from warm, photo- 
or shock-ionized gas at ~ 10000 K (Kobulnicky & Skill- 
man 1997). 
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However, although the metal abundances in the hot 
regions are quite large, reaching also extrasolar values 
(cf. Fig. 14), in our models the majority of the metals 
is in the cold gas phase. Thus the abundances in Ta- 
ble 4 are essentially the same as those in the cold gas, 
and their comparison with the observed metallicities is 
meaningfull. Note that, in making such a comparison, 
we suggest that the present time burst in IZwl8 can be 
responsible of the observed chemical enrichment in the 
H II regions. This is in agreement with the arguments 
illustrated in section 3.2, and at variance with previ- 
ous suggestions of different authors (see e.g. Larsen et 
al. 2000). Recent observations (Pettini & Lipman 1995; 
Van Zee et al. 1998), although uncertain, indicate an 
oxygen abundance in the Hi regions of IZwl8 compara- 
ble with that in the H II regions, in agreement with our 
predictions. 

Table 4 reports the abundances of our models after 
31 Myr. The lifetime of a 8 Mq star is approximatively 
29 Myr, according to eq. (14). Therefore, since at Z = 
secondary N is not produced and primary N from mas- 
sive stars is negligible, only for ages larger than 29 Myr 
we can expect some N which is the one produced in 
a primary fashion by IMS during the third dredge-up 
episode. As one can see from Table 4, the abundances 
and abundance ratios predicted by the yields for Z = 
and cxrv = 1-5 (model M1B) are in good agreement 
with those measured in IZwl8, thus we could conclude 
that the abundances in this galaxy are compatible with 
only one burst, the first, but only if the burst age is 
of the order of 31 Myr. In fact, for times shorter than 
that the N abundance is too low and for times larger 
the agreement worsens. 

The abundance of C and particularly the predicted 
C/O ratio in model M1B is in very good agreement with 
observations at variance with previous works (Kunth et 
al. 1995). The difference between the low abundance 
of 12 C predicted by Kunth et al. (1995) and here is, 
in our opinion, due to the fact that the total amount 
of stars produced there was smaller (M bu rst = 2 x 10 5 
Mq) and less in agreement with the observations than 
that produced here (Mburst = 6 x 10 5 M ), therefore 
we predict higher abundances for all the elements. In 
addition, the yields for massive stars used here are dif- 
ferent from those used in Kunth et al. (1994) (those of 
Woosley 1987). It is also worth mentioning that the es- 
timated ages for the present burst in IZwl8 are between 
15 and 27 Myr (Martin 1996) in good agreement with 
our suggestion, although other authors suggest ages as 
short as 5 Myr (Izotov et al. 1999; Stasihska & Schaerer 
1999). 

In order to see if we can exclude a previous burst 
besides the present one in IZwl8, as suggested by pre- 
vious papers (Aloisi et al. 1999; Kunth et al. 1995), or 
a recent burst coupled with a low but continuous star 
formation (Legrand 1999), we computed the expected 
ISM abundances for a preenriched gas with Z = 0.01 
Zq and we show in Table 4 the abundances for this case 
at an age of 31 Myr (cases C and D). The results for the 
C case show that at an age of 31 Myr the abundance 
of oxygen is too high; the results at the end of the sim- 
ulation (375 Myr for model Ml) give better values for 
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Figure 12. Abundance gradients for oxygen (solid line) and 
iron (dotted line) after 31 Myr, for the model M1B. 

oxygen and N/O but they predict a too high C/O. If one 
assumes then Z = 0.01 Z Q and any = 1.5 (case D), the 
agreement worsens at any age since one predicts a too 
high N/O ratio while the rest is practically unchanged. 
Therefore, we have two considerations: first, the single 
first burst hypothesis seems to give the best agreement 
with observations as long as primary N production in 
IMS is considered; second, we cannot really exclude a 
previous burst before the present one, or at least we 
cannot exclude a burst which enriched only sligthly the 
ISM. In other words, the preenrichment should be less 
than Z — 0.01 Zq. From the previous discussion, it 
arises that the best age for the burst in IZwl8 should 
be around 31 Myr. 

At this age the abundance gradient in a region of 
600 pc is almost flat, at least if the galaxy is seen edge- 
on (see Fig. 12), in agreement with what is observed 
in IZwl8 (Legrand 1999). Actually, if a bipolar-shaped 
expanding bubble is present in IZwl8, the inclination 
of the symmetry axis with respect to the normal to the 
observer would be very small (Martin 1996, suggested 
an inclination of 10°). 

Finally, in Fig. 13 we show the predicted [q/Fc] 
vs. time and vs. [Fe/H] for the gas inside and outside 
the galaxy, corresponding to the results of Models M1B, 
M2B and M3B. The interesting feature of this figure is 
that the [a/Fe] ratios in the gas outside the galaxy are 
lower than those in the gas inside the galaxy. This is 
due to the fact that Fe, in particular that produced by 
type la SNe, is lost more efficiently than a elements. 

We shall also note in these two figures the peculiar 
behaviour of silicon relative to the other a-elements. 
Silicon is, in fact, synthetized by SNela at variance with 
what happens for O and Mg, which are mainly produced 
by SNell (see Gibson et al. 1997). This is the reason for 
the relatively low and flat [Si/Fe] ratio as a function of 
[Fe/H] . 



5 CONCLUSIONS 

We have studied the dynamical and chemical evolution 
of a dwarf galaxy as due to the effect of a single, instan- 
taneous, point-like starburst occurring in its centre. We 
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Figure 13. Predicted [a/Fe] vs. time and vs. [Fe/H] for models M1B, M2B, M3B (first, second and third raw respectively) 
for both expelled gas and ISM 



adopted galactic structural parameters which resemble 
those of IZwl8, the most unevolved dwarf blue compact 
galaxy known locally. We ran different models, which 
differ for the burst luminosity and the fSM mass. We 
considered the mass and energy inputs from the single 
low and IMS, the SNell and the SNela (white dwarfs 
in binary systems) and we followed the evolution of the 
gas and its chemical abundances (H, He, C, N, O, Mg, 
Si and Fe) in space and time for several hundreds Myr 
from the burst. 

Our results can be summarized as follows: 



(i) The starburst can inject enough energy into the 
ISM to trigger a metal enriched galactic wind: the met- 
als synthetized and ejected through supernova explo- 
sions leave the galaxy more easily than the unprocessed 
gas. This result is not new since it has already been 
suggested by previous works (e. g. MF and DB). How- 
ever, our new result is that the SNIa products have the 
largest ejection efficiency (more than the products of 
type II SNe), with the consequence that the [a/Fe] ra- 
tios in the gas outside the galaxy are predicted to be 
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Figure 14. Abundance evolution in hot and cold regions for the model M1B. The threshold temperature is 2 X 10 4 K. 



lower than those inside. This is due to the fact that 
SNela produce a substantial fraction of iron. 

(ii) The energy injection in the ISM by SNII has a 
rather low efficiency. Instead, the energetic contribution 
of SNela, in spite of their relatively small number (a to- 
tal of 240 SNela against 4800 SNell), has important 
consequences in the dynamical behaviour of the galaxy. 
Since the SNIa explosions occur in a medium heated and 
diluted by the previous activity of SNell, the thermal 
energy of the explosions is easily converted into kinetic 
energy and so the gas reaches quickly the outer regions 
of the bubble along the galactic chimney. DB showed 
that after the end of SNII activity, a fraction of the 
gas tends to recollapse toward the central region of the 
galaxy, achieving the threshold density for a new star 
formation event (iV(Hl) > 10 21 cm~ 2 , Skillman et al. 
1988; Saito et al. 1992) in ~ 0.5-1 Gyr. With the ener- 
getic contribution of SNela, it is not possible to reach, 



at least for the time considered in our simulation, such 
a threshold in column density (it can be obtained only 
after many Gyr). 

(iii) One single burst, occurring in a primordial gas, 
with an age of ~ 31 Myr reproduces quite well both 
the dynamical structures and the abundances in IZwl8. 
This value is consistent with other independent age es- 
timates for the burst (Martin 1996). From the nucle- 
osynthetic point of view the age of 31 Myr ensures that 
there is enough time for the primary N from IMS to be 
produced and ejected. The adopted yields for massive 
stars (Woosley & Weaver 1995) include some primary 
N, but its amount is negligible. This result suggests that 
IZwl8 is probably experiencing its first major burst of 
star formation, although we cannot exclude a previous 
burst (see Aloisi et al. 1999) of moderate intensity, which 
enriched the gas to a metallicity Z < 0.01 Zq. 

(iv) At variance with previous studies we find that 
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the majority of metals (in mass) are found in the cold 
gas. In fact, mainly because of the low SNII efficiency, 
the wind superbubble remains several hundreds of Myr 
inside the galaxy before the break out occurs. Moreover, 
the superbubble becomes radiative after a few Myr, and 
most of the SNII ejecta cool without leaving the galaxy 
(SNIa ejecta, instead, are shown to be vented more eas- 
ily). Given the relatively short mixing time, the abun- 
dances predicted by our models for the cold gas are those 
that should be compared with the abundances observed 
in IZwl8. Actually, they are in very good agreement 
with the observed ones. This result supports the com- 
mon assumption made in chemical evolution models of 
an instantaneous mixing. 

Future improvements of this work will include mod- 
els with continuous bursts, sequential bursts, as well as a 
more detailed study of the formation of the H II regions 
(Recchi et al. in preparation). 



ACKNOWLEDGEMENTS 

We are grateful to Guillermo Tenorio-Tagle and Andrea 
Ferrara for useful suggestions and discussions. We also 
thank the referee whose suggestions improved the paper. 
We acknowledge financial support from the Italian Min- 
istry for University and for Scientific and Technological 
Research (MURST). 



REFERENCES 

Aloisi A., Tosi M., Greggio L., 1999, AJ, 118, 302 
Anders E., Grevesse N., 1989, Geochim. Cosmochim. Acta, 
53, 197 

Bedogni R., D'Ercole A., 1986, A&A, 157, 101 
Begelman M.C., Fabian A., 1990, MNRAS, 244, 26p 
Begelman M.C., McKee C.F., 1990, ApJ, 358, 375 
Borkowski K.J., Balbus S.A., Fristrom C.C., 1990, ApJ, 355, 
501 

Bradamante F., Matteucci F., D'Ercole A., 1998, A&A, 337, 
338 

Breitschwerdt P., Kahn F.D., 1988, MNRAS, 235, 1011 
Ciotti L., D'Ercole A., Pellegrini S., Renzini A. 1991, ApJ, 
376, 380 

Cowie L.L., McKee C.F., 1977, ApJ, 215, 213 
Cowie L.L., McKee C.F., Ostriker J. P., 1981, ApJ, 247, 908 
Dekel A., Silk J., 1986, ApJ, 303, 39 
D'Ercole A., Brighenti F., 1999, MNRAS, 309, 941 (DB) 
Dcvost D., Roy J.R., Drissen R., 1997, ApJ, 482, 765 
De Young D.S., Gallagher J.S., 1990, ApJ, 356, L15 
De Young D.S., Heckman T.M., 1994, ApJ, 431, 598 
Dufour R.J., Esteban C, Castaneda H.O., 1996, ApJ, 471, 
L87 

Dufour R.J., Hester J. J., 1990, ApJ, 350, 149 

Garnett D.R., Dufour R.J., Peimbert M., Torres-Peimbert 
S., Shields G.A., Skillman E.D., Terlevich E., Terlevich 
R.J., 1995, ApJ, 449, L77 

Garnett D.R., Skillman E.D., Dufour R.J., Shields G.A., 
1997, ApJ, 481, 174 

Gibson B.K., Loewenstein M., Mushotzky R.F., 1997, MN- 
RAS, 290, 623 

Gibson B.K., Matteucci F., 1997, ApJ, 475, 47 

Greggio L., Renzini A., 1983, A&A, 217, 222 

Hunter D.A., Thronson H.A., 1995, ApJ, 452, 238 



Israel F. P., 1988, A&A, 194, 24 

Izotov Y.I., Paparedos P., T huan T.X.. Fricke K.J.. Foltz 

G.B., Guseva N.G., 1999, ^stro-ph/9907082 
Kahn F.D., Breitschwerdt P., 1989, MNRAS, 242, 209 
Kobulnicky H.A., Skillman E.D., 1997, ApJ, 489, 636 
Koo B.C., McKee C.F., 1992, ApJ, 388, 9 
Kunth D., Lequeux J., Sargent W.L.W., Viallefond F., 1994, 

A&A, 282, 709 
Kunth D., Matteucci F., Marconi G., 1995, A&A, 297, 634 
Lar scn T.I.. Somme r-Larsen J., Pagel B.E.J., 2000, 

|astro-ph/0005249| 
Larson R.B., 1974, MNRAS, 169, 229 
Legrand F., 1999, Ph. D. Thesis, Univ. Paris 6 
Leitherer C, Schaerer D., Goldader J.D., Gonzales Delgado 

R.M., Robert C, Kune R.F., de Mello D.F., Devost D., 

Heckman T.M., 1999, ApJS, 123, 3 
Lequeux J., Viallefond F., 1980. AfcA. 91. 269 
Lin D.N.C., Murray S.D., 2000, nstro-ph/0004055 
Loewenstein M., Mathews W.G., 1987, ApJ, 319, 471 
Mac Low M.-M., Ferrara A., 1999, ApJ, 513, 142 (MF) 
Mac Low M.-M., McCray R., 1988, ApJ, 324, 776 
Mac Low M.-M., McCray R., Norman M.L., 1989, ApJ, 337, 

141 

Marconi G., Matteucci F., Tosi M., 1994, MNRAS, 217, 391 
Marlowe A.T., Heckman T.M., Wyse R.F.G., Schommer R., 

1995 ApJ, 438, 563 
Martin C.L., 1996, ApJ, 465, 680 

Mas-Hesse J.M., Kunth D., 1996, in The Interplay between 
Massive Star Formation, the ISM and Galaxy Evolution, 
ed. D. Kunth, B. Guideroni, M. Heydari-Malayeri & T.X. 
Thuan (Gif-Sur-Yvette: Ed. Frontieres), 401 
Mathews W.G., Bregman J.N., 1978, ApJ, 224, 308 
Matteucci F., 1986, MNRAS, 221, 911 
Matteucci, F., 1996, Fond. Cosm. Phys., 17, 283 
Matteucci F., Chiosi C, 1983, A&A, 123, 121 
Matteucci F., Tosi M., 1985, MNRAS, 217, 391 
McKee C.F., Begelman M.C., 1990, ApJ, 358, 392 
McKee C.F., Van Buren D., Lazareff B., 1984, ApJ, 278, 
L115 

Meurer G.R., 1991, Proc. Astron. Soc. Australia, 9, 98 
Meurer G.R., Freeman K.C., Dopita M.A., Cacciari C, 1992, 
AJ, 103, 60 

Nomoto K., Thielemann F.K., Yokoi K., 1984, ApJ, 286, 644 
(NTY) 

Olofsson K., 1995, A&A, 293, 652 
Padovani P., Matteucci F., 1993, ApJ, 416, 26 
Pettini M., Lipman K., 1995, A&A, 297, L63 
Pilyugin L.S., 1992, A&A, 260, 58 
Pilyugin L.S., 1993, A&A, 277, 42 

Portinari L., Chiosi C, Bressan A., 1998, A&A, 334, 505 
Reimers D.: 1975, Mem. R. Sci. Liege 6 erne Ser., 8, 369 
Renzini A., Voli M., 1981, A&A, 94, 175 (RV) 
Rieschick A., Hensler G., 2000, in Cosmic Evolution and 

Galaxy Formation: Structure, Interactions and Feedback, 

ASP Conferences Series, eds. J. Franco, E. Terlevich, O. 

Lopez Cruz and I. Aretxaga 
Roy J.-R., Kunth D., 1995, A&A, 294, 432 
Sait5 M., Sasaki M., Ohta K., Yamada T., 1992, PASJ, 44, 

593 

Salpeter E.E., 1955, ApJ, 121, 161 
Sandage A., Binggeli B., 1984, AJ, 89, 919 
Silich S.A., Tenorio-Tagle G., 1998, MNRAS, 299, 249 
Skillman E.D., 1997, Rev. Mex. SC, 6, 36 
Skillman E.D., Bender R., 1995, Rev. Mex. SC, 3, 25 
Skillman E.D., Kennicut R.C.J., 1993, ApJ, 411, 655 
Skillman E.D., Kennicut R.C.J., Hodge P.W., 1989, ApJ, 
347, 875 



© 0000 RAS, MNRAS 000, 000-000 



Dynamical and chemical evolution of gas-rich dwarf galaxies 

Skillman E.D., Terlevich R., Teuben P.J., van Wocrdcn H., 

1988, A&A, 198, 33 
Slavin J.D., Shull J.M., Begelman M.C., 1993, ApJ, 407, 83 
Smith M.D., Smarr L., Norman M.L., Wilson J.R., 1983, 

ApJ, 264, 432 

Spitzcr L., 1956, Physics of Fully Ionized Gases (New York: 

Interscicncc) 
Stasihka G., Schaerer D., 1999, A&A, 351, 72 
Stone J.M., Norman M.L., 1992, ApJS, 80, 573 
Stone J.M., Norman M.L., 1993, ApJ, 413, 198 
Strickland D.K., Stevens L.R., 1999, MNRAS, 306, 43 
Suchkov A.A., Balsara T.S., Heckman T.M., Leiterher C, 

1994, ApJ, 430, 511 
Tenorio-Taglc G., 1996, AJ, 111, 1641 
Vader J.P., 1986, ApJ, 379, LI 

Van Zee L., Westpfahl D., Haynes M.P., Salzer J. J., 1998, 
AJ, 115, 1000 

Viallcfond F., Lequcux J., Comte G., 1987, in Starburst and 
Galaxy Evolution, ed. T.X. Thuan, T. Montmerle & J. 
Tran Thanh Van (Gif-sur-Yvette: Ed. Frontieres), 139 

Weaver R., McCray R., Castor J., Shapiro P., Moore R., 
1977, ApJ, 218, 377 

Woosley S.E., 1987, in Nucleosynthesis and Chemical evolu- 
tion, 16th Saas-Fec Advanced Course, Geneva Observa- 
tory, 1 

Woosley S.E., Weaver, T.A., 1995, ApJS, 101, 181 (WW) 



© 0000 RAS, MNRAS 000, 000-000 



This figure "recchi5a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0002370v2 



This figure "recchi5b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0002370v2 



This figure "recchi6a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0002370v2 



This figure "recchi6b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0002370v2 



This figure "recchi7a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0002370v2 



This figure "recchi7b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0002370v2 



This figure "recchi8a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0002370v2 



This figure "recchi8b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0002370v2 



